基于自适应网格的钝头体激波诱导燃烧现象数值模拟研究

钝头体激波诱导燃烧是爆震研究的一个基本问题。针对化学恰当量比的H2/Air预混气在Ma=4.79和Ma=6.46时的激波诱导燃烧现象开展数值模拟研究，采用基于有限体积法的块结构自适应网格加密程序AMROC对带化学反应源项的轴对称Euler方程解耦求解，考察了数值模拟中不同形式的MUSCL重构格式、限制器类型以及化学反应机理等重要因素对模拟结果的影响。结果表明，程序能够根据设定的加密判据较好地实现网格自适应加密，减小总网格量，实现高效数值模拟。通过与实验数据的对比，表明非定常激波诱导燃烧算例的准确程度不仅取决于化学反应机理，也取决于限制器类型，而采用两种不同形式的MUSCL重构格式获得的振荡频率则几乎一致，与试验结果的误差分别为1.17%和0.97%。模拟对比经典的Jachimowski机理和近年来新发展的几种包含压力相关反应步的氢/氧反应机理，模拟结果表明：对于Ma=4.79时的非定常激波诱导燃烧模拟，经典的Jachimowski机理仍然是能够给出与实验结果最接近的反应机理；而对于Ma=6.46时的定常激波诱导燃烧模拟，几种反应机理均能给出与实验吻合较好的结果。     

碳纳米管阵列推力器放电特性的数值模拟研究

基于直接蒙特卡洛(DSMC)算法、浸入式有限元算法(IFE)、粒子云及蒙特卡洛碰撞(PIC-MCC)算法以及Ammosov-Delone-Krainov(ADK)隧穿电离模型建立了碳纳米管阵列推力器工质气体和CNT场电离过程的三维仿真模型。利用文献中的实验数据对本文仿真模型的有效性进行了验证,结果表明本文仿真模型与实验结果基本一致。针对碳纳米管阵列推力器的放电特性进行了研究,并重点分析了离子与原子碰撞过程对碳纳米管阵列推力器放电特性的影响。结果表明CNT可以极大地增强局部电场,稳态时场增强因子约为1308;与原子的碰撞使得离子在CNT附近更加集中,降低了局部电场,进而降低了其场增强能力。     

不同直径圆球诱导燃烧的振荡机制与频率特性

为了分析不同直径圆球诱导振荡燃烧的规律，并揭示圆球大小在振荡燃烧现象中所发挥的深层次作用，本文采用二维轴对称Euler方程和基元反应模型，对不同直径的圆球在H2/air预混气体中诱导振荡燃烧的现象开展数值模拟研究。研究发现，振荡频率并不是简单地随直径增大而逐渐从高频向低频连续过渡，而是存在两次突变，形成了超高频、高频以及低频三种振荡燃烧模态。在两种模态间过渡时，振荡达到稳定状态前，会存在一段双频耦合的振荡阶段。三种不同振荡燃烧模态的产生是受到了不同振荡机制的作用，而两种模态间过渡时的双频耦合现象则是两种机制相互竞争的结果。     

偏心撞击对撞击式喷嘴雾化特性的影响

为研究偏心撞击对撞击式喷嘴雾化特性的影响,建立了求解自燃推进剂冷态射流撞击雾化过程的数值模拟方案,计算了不同偏心度条件下的射流撞击雾化过程。采用树形自适应加密算法直接求解不可压Navier-Stokes方程组,由分段线性的流体体积(VOF)方法对流体界面进行捕捉。结果表明偏心撞击会导致雾场发生偏转,当无量纲偏心度E为1/8时,雾场偏转角度约为9.2°,应控制加工偏差小于该值。随着偏心度的增大,液膜的偏转角度增大,理论推导得到的液膜偏转角度要小于数值计算得到的液膜偏转角度。正心撞击时燃料与氧化剂流强峰值接近,雾场的流强分布呈单峰分布。当发生偏心撞击时,由于燃料与氧化剂部分射流未参与撞击导致流强峰值出现交错,雾场的流强分布呈双峰分布,混合比的空间分布发生较大改变。正心撞击时撞击点下游液滴的速度分布近似呈轴对称分布,而偏心撞击之后的速度分布则呈中心对称分布。偏心撞击导致的射流动量损失使得雾化性能变差,当无量纲偏心度E为1/8时,一甲基肼(MMH)的Sauter平均直径增大约4.8%,四氧化二氮(NTO)的Sauter平均直径增大约5.8%。      

跨介质航行器波浪环境入水流场演变和运动特性研究

跨介质航行器是一种既可以在空中飞行又可以在水下潜航的新概念航行器,基于仿生学原理,提出一种通过改变外形实现水空介质跨越的航行器模型,通过入水试验装置和计算流体动力学方法,对航行器带攻角从空气到水的介质跨越过程进行了试验和数值仿真研究,得到了跨介质入水过程航行器的运动姿态和入水空泡形态,并通过数值仿真得到了航行器的升力、阻力、速度和加速度演化规律。同时基于数值模拟方法对有波浪情况和静水情况下航行器入水过程空泡演变以及运动特性进行对比。结果表明:提出的航行器构型在水中具有较好的姿态调整能力,波浪的有无和波高的不同都会对航行体入水运动特性造成影响。     

Initiation characteristics of oblique detonation waves from a finite wedge under argon dilution

In this paper, the flow field characteristics of Oblique Detonation Waves (ODWs) induced by a finite wedge under argon dilution are studied by solving the Euler equations with a detailed chemical model of hydrogen and air. First, the effects of the expansion waves, argon concentration, geometric parameters, and Mach number on the ODW are discussed. The results show that the changes of these parameters may make the oblique detonation not be initiated. Then, the ODW initiation criterion of the finite wedge is summarized, as the characteristic length of the induction zone L_C and the characteristic length of the oblique wedge L_W meet the condition L_C/L_W < 1, the initiation of the ODW occurs; otherwise, it does not occur. What’s more, the Constant Volume Combustion (CVC) theory is applied to study the characteristic length of induction zone. It is found that CVC theory is more suitable for the “smooth transition” type of ODW flow field, the theoretical and numerical characteristic length in induction regions are in good agreement. This work is of great significance for the design of oblique detonation engines for hypersonic vehicles.     

燃烧室头部激励的等离子体强化燃烧特性实验研究

等离子体助燃是一种新型的强化燃烧技术，近年来受到国内外学者的广泛关注。本文开创性地研制了基于旋转滑动弧等离子体的强化燃烧头部，建立了某型航空发动机三头部燃烧室实验件的等离子体助燃实验平台，验证了该等离子体强化燃烧技术应用于型号发动机燃烧室的可行性。实验研究等离子体助燃在不同余气系数和不同输入电压条件下对平均出口温度、燃烧效率、温度分布系数以及熄火边界的影响。实验结果表明，与正常燃烧相比，施加等离子体助燃后的燃烧效率有明显的提高，在输入电压为U0=240V，余气系数为 α=0.8时，等离子体助燃的燃烧效率提高3.24%。实施等离子体助燃后，燃烧室出口温度分布场分布得到明显的改善，在富油工况α=0.8，出口温度分布系数减少39.8%。等离子体助燃输入电压越高熄火边界扩展程度越明显，相比于正常工况条件下，等离子体助燃U0=240V的熄火边界扩宽了7.34%。     

直升机新型旋翼翼型气动特性与布局分析

利用结构网格计算流体力学(Computational fluid dynamics ,CFD)的翼型气动特性分析方法开展旋翼翼型气动特性计算。通过RAE282，NACA0012，OA212，OA207等翼型压力分布、升力和阻力等特性计算结果与试验结果的对比分析，验证了计算方法的准确性，并进一步完成了HF系列旋翼翼型的气动特性计算。基于翼型的气动特性，采用时间步进自由尾迹的旋翼气动性能分析方法开展旋翼桨叶翼型的气动布局优化设计，对悬停和前飞条件下的旋翼开展计算分析，得到两种条件下的旋翼气动特性。而后通过本文建立的优化方法开展旋翼翼型布局优化设计。     

非恒定增压式微尺度拉瓦尔喷管流动特性

针对微尺度喷流在航天器运动状态切换时出现的非恒定增压变化,采用直接模拟蒙特卡洛(DSMC)方法对阶跃式增压和线性式增压两种模式下的微尺度拉瓦尔喷管流场进行了模拟,并对其变化过程中的流动特性进行了对比分析。结果显示：阶跃式增压会导致流动特性出现较大幅值的峰谷式波动,而线性式增压下的流动特性则呈现出线性变化的特点;黏性力对微尺度喷流的非恒定增压变化产生了重要的黏滞作用,在喉部扩张段至出口的流场中尤为明显;在设定的条件下,阶跃式增压过程中喷流产生的总冲量较线性式增压高59.5%,质量流量高74.7%,单位工质提供的冲量低约8.6%,波动性也高于线性式模型,阶跃式增压适用于系统需要较大推力改变运动状态且推进剂充足的情况,而线性式增压在系统精确微调或需要推进剂产生更高效能时具有明显的优势。     

针栓式推力室冷却特性试验研究

针栓式发动机具有推力调节简单、声学燃烧稳定性好和成本低廉等一系列优点,但目前对针栓式推力室的冷却特性掌握还不够深入,研制过程中曾多次出现因推力室冷却组织不好而烧蚀的问题。采用针栓式推力室试验件,对身部的冷却特性进行了深入研究,结果表明:推力室圆柱段前端的温度较低,无需采用专门的冷却措施;末端内侧气壁温度会达到约1650℃,局部存在明显的烧蚀现象,必须采取有效的热防护措施;下游冷却液膜量的变化对圆柱段的冷却特性影响较小。